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In compressed supersymmetry, a light top squark naturally mediates efficient neu- 

tralino pair annihilation to govern the thermal relic abundance of dark matter. I 

study the LHC signal of same-sign leptonic top-quark decays from gluino and squark 

production, which follows from gluino decays to top plus stop followed by the stop 

' r~|'. decaying to a charm quark and the LSP in these models. Measurements of the num- 

*^\ bers of jets with heavy-flavor tags in the same-sign lepton events can be used to 

C^l confirm the origin of the signal. Summed transverse momentum observables provide 

r~| \ an estimate of an effective superpartner mass, which is correlated with the gluino 

mass. Measurements of invariant mass endpoints from the visible products of gluino 

CN ■ decays do not allow direct determination of superpartner masses, but can place con- 

JI^ ■ straints on them, including lower bounds on the gluino mass as a function of the 

CN . top-squark mass. 
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I. INTRODUCTION 

If supersymmetry [li] is the solution to the hierarchy problem associated with the small ratio of 
the electroweak scale to the Planck scale, then some of the superpartners should be discovered at 
the impending CERN Large Hadron Collider (LHC). While the essential idea of supersymmetry 
as a symmetry connecting fermion and boson degrees of freedom is quite predictive, the unknown 
features of the supersymmetry breaking mechanism allow for a diverse variety of possibilities for 
the LHC signals of superpartner production and decay [2, y] • 

The purpose of this paper is to study some of the distinctive LHC signals particular to the 
"compressed supersymmetry" scenario proposed in ref. [4], which is motivated both by the su- 
persymmetric little hierarchy problem and by the cold dark matter relic abundance obtained by 
WMAP, SDSS and other experiments [5|, l6(. This model scenario follows from assuming that the 
ratio of the running gluino and wino mass parameters, M3/M2, is smaller than 1 near the GUT 
scale, unlike the assumption of the well-studied "minimal supergravity" (mSUGRA) framework. 
(Models with non- unified gaugino masses have recently attracted renewed interest, see for example 
m- 23], due to their ability to incorporate novel LHC phenomenology and dark-matter physics.) As 



a result, the ratio of the physical masses of the heaviest and the lightest superpartners is much less 
than in mSUGRA, because the gluino mass feeds into the other superpartner masses by renormal- 
ization group evolution. Another characteristic feature is that the pair annihilation of the lightest 
supersymmetric particles (LSPs) in the early universe can naturally proceed dominantly through 
NiNi — > ti, mediated by top-squark exchange. This is due to the fact that the stop-LSP mass 
difference is naturally not too large, particularly in models that have enough top-squark mixing to 
evade the Higgs scalar boson mass bound from LEP. 

In compressed supersymmetry, the superpartner masses are typically all less than 1 TeV in the 
available parameter space. This means that the initial evidence for supersymmetry should follow 
quickly from the classic jets with missing transverse energy (i?™'^^) signal, as soon as the systematic 
difficulties associated with understanding missing energy as manifested in the LHC detectors are 
conquered. One can then turn attention to those features of the signal that might distinguish it 
from the usual mSUGRA models. In much of the parameter space in compressed supersymmetry 
that predicts the observed thermal relic abundance of dark matter, the mass difference between 
the top squark and the neutralino LSP A^i is less than 85 GeV, so that the flavor-preserving decays 
including ii — > tNi and ti — > bCi and ti — > bWNi are kinematically forbidden. In this paper, I 
will assume that the flavor-violating 2-body decay 

h ^ cNi (1.1) 



dominates over the remaining possibility 3J, |35|], the 4-body decay ti -^ bff'Ni. (For more on 



this assumption, see the next section.) Assuming the top squark is lighter than the gluino, one 
has: 

(til (50%) , , 

g^l ! (1.2) 

[ tti (50%), 



due to the Majorana nature of the gluino, leading to 



PP^ 99 



ttccNiNi (25%) 

ttccNiNi (25%) (1-3) 

ttccNiNi (50%). 



When both of the same-sign top quarks of the first two cases decay leptonically, one obtains a 
distinctive detector signal of two same-sign leptons, two potentially 6-tagged jets, two or more 
additional jets, and missing energy from the LSPs and neutrinos: 



pp^i^eHbjj + Ef'^\ (1.4) 

This is a special case of the well-known same-sign dilepton signature for Majorana gluino (or 
gaugino) production in supersymmetry [3g]. This LHC signal for gluino pairs was proposed and 
studied in some detail, in the context of models with much lighter top squarks (m? < mt), in 



ref. 37|. Adding to this signal in compressed supersymmetry will be events in which squarks are 
produced, giving extra jets in the final state when they decay to the gluino. The presence of same- 
sign leptons provides for a strongly suppressed Standard Model background compared to other 
missing energy signals, and this is further aided by requiring two 6-tagged jets. 

In this paper, I will consider the properties of the LHC events that conform to this signal 
in compressed supersymmetry. (Other studies of the LHC phenomenology of compressed super- 



symmetry are found in [38|- 391].) Section [III defines a model line for study, a one-parameter slice 
of model space with the free parameter corresponding to the gaugino mass scale. I also discuss 
some of the prominent properties of the superpartner mass spectrum of this model line that make 
it qualitatively different from mSUGRA models. Section IIIII describes an event selection for the 
i^f^bbjj + i?™^^, and the features of the resulting signal events. Section IIVI considers mass- 
estimating observables based on the scalar sum of transverse momentum of detector objects, while 
section |V] studies kinematic endpoints of the invariant masses of visible products of the gluino 
decay. Section IVTl contains some concluding remarks. 



II. A COMPRESSED SUPERSYMMETRY MODEL LINE 

One simple realization of compressed supersymmetry is obtained by supposing that the running 
bino, wino, and gluino masses are parameterized at Mgut by: 

Ml = mi/2(l + C24), (2.1) 

M2 = mi/2(l + 3C24), (2.2) 

Ms = mi/2(l - 2C24), (2.3) 

corresponding to an F-term source for supers ym metry breaking in a linear combination of the 



singlet and adjoint representations of SU{5) [7[- 10(]. Merely for simplicity, I also assume a common 



scalar mass ttt-q and scalar trilinear coupling Aq, both at Mqut- The other parameters defining 



the model are tan/5 and the phase of the /i parameter, which is taken to be real. I use SDFTSUSY 
2.0.11 [4(J] to generate the superpartner spectrum. To define the model line for study here, let: 

C24 = 0.21, Ao/Mi = -1, tan/3 = 10, /x > 0, (2.4) 

with Ml (or equivalently mx/2) taken as the single varying parameter of the model line. (Here and 
in the following, Mi is used to denote the running bino mass parameter at Mqut ) not at the elec- 
troweak scale.) For each value of Mi, the parameter niQ is obtained by imposing as a requirement 
that the predicted dark matter relic abundance (obtained using the program micrOMEGAs 2.0.1 
411 1) satisfies Odm^^ = 0.11 [5|, la]- The resulting values of rriQ are not too large, ranging from 210 



to 380 GeV (and always less than the wino and bino masses at the GUT scale) for the model line 
when the physical gluino mass is less than 1 TeV. 

With strictly flavor-conserving boundary conditions for the soft super symmetry-breaking inter- 
actions at the GUT scale, the 2-body decay ti — > cNi and the 4-body decay ti -^ hf fN-\ would 



have roughly comparable partial decay widths on this model line. Using SDECAY 4^], one finds 
that BR(ii — > bff'Ni) would range from a few percent (for a small mass difference m^ — w-jy ~ 30 
GeV), to nearly 90% (for the largest mass difference of about 70 GeV).^^ However, the strict minimal 
flavor violation assumption on which this is based is notoriously unmotivated by theory, except in 
models with special features like gauge-mediated supersymmetry breaking. A small amount of non- 
minimal flavor violation results in the 2-body decay ti — > cA'^i dominating, as assumed here. Writing 
the effective stop-charm-neutralino interaction at the weak scale as £ = —t\(jjLCLNi + UrcrNi), 
I find that even for the worst-case point on the model line, BR(ti — > cA'^i) > 95% provided that 
(yi + yW > 8 X 10"'^. This would follow, for example, from a small off-diagonal right-handed 
up squark squared mass parameter at the GUT scale m?, - /rrig > 0.007, with no danger of conflict 
with present flavor experiments. 

The superpartner and Higgs boson mass spectrum for a representative point on the model 
line with Mi = 500 GeV is shown in flgure [TJ The ratio of masses of the heaviest and lightest 
superpartners in this model is 3.6, almost a factor of 2 smaller than is obtainable in mSUGRA 
models even with small rriQ.The neutralino LSP is heavier than the top quark, allowing NiNi -^ ti. 
The top squark is the next-to-lightest superpartner. Another distinctive feature is that Ci and 
A''2,3 are higgsino-like states, due to the fact that fi is only 361 GeV, again much smaller than 
found in mSUGRA models for comparable gluino and squark masses. (The relatively small value 
of |/x| is a sign of the reduced flne-tuning found in models with a small ratio M3/M2 at the GUT 
scale as pointed out long ago in jll|].) The heaviest neutralino and chargino states N4 and C2 are 
wino-like. The sleptons turn out to be too heavy to play a signiflcant role in LHC physics. 

The gluino always decays to top and stop in this model, and the stop always decays to a charm 
quark and LSP. Other important decay modes are, for left-handed squarks: 

UL^ug {71%), dC2 (13%), uiV4 (6%), dCi {6%), (2.5) 

di^dg (73%), nC'2 (14%), dN^ (7%), uCi (3%), (2.6) 



^ The most important contribution to the 4-body decay partial width in this model line comes from the Feynman 
graph with virtual W and top-quark exchange, unlike in most other models considered in the literature (see, for 
example, [35l. |43|. |44|]). where diagrams with W and chargino or slepton and chargino exchange dominate. 



FIG. 1: The mass spectrum for a sample 
point on the model line described in the 
text, with Ml = 500 GeV at the GUT 
scale, C24 = 0.21, Aq/Mx = -1, mo = 
314 GeV, ^Ji = 361 GeV, and tan/3 = 
10. The columns contain, from left to 
right, Higgs scalar bosons, ncutralinos, 
charginos, the gluino, first and second 
family squarks and sleptons, and third 
family squarks and sleptons. 
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and for right-handed squarks: 



UR^uNi (92%), ug (5%), UN2 (3%), 
dR^dNi (85%), dg (12%), diVs (3%). 



(2.7) 
(2.8) 



Thus left-handed squarks are a plenteous source of gluinos, while right-handed squarks mostly 
decay directly to the LSP. Subdominant decays produce some neutralinos and charginos, which 
nearly always decay into on-shell W, Z, and h bosons or through top squarks. For the higgsino-like 
states: 



N2 -^ hNi (90%), ZNi (10%), 
iVg ^ hNi (97%), ZNi (3%), 
Ci ^ bii (91%), WNi (9%), 



(2.9) 
(2.10) 
(2.11) 



and for the heavier, wino-like, states: 



N4 -^ WCi (51%), 
C2 -^ WN2 (30%), 



hN2 (20%), ZN3 (20%), tti 

WN3 (21%), ZCi (25%), hCi (21%) 



(2.12) 
(2.13) 



An important consequence of these decays is that one cannot find dilepton mass edges of the type 
used in 45|]- 49|] to obtain information about the superpartner mass spectrum. The only isolated 



leptons come from on-shell W and Z decays, since two-body spoiler decays are always allowed. 
Furthermore, sleptons completely decouple from the cascade decays, because they are too heavy. 
These features are qualitatively maintained along the entire model line. 

Varying Mi, one finds that Mi > 417 GeV at the GUT scale is required to satisfy the LEP 
bound on the Higgs mass, taken here to be m/j > 113 GeV because of the theoretical uncertainty 
on the Higgs mass prediction. As Mi increases, the gluino and LSP masses increase approximately 
in direct proportion, while the top-squark mass stays between 30 and 70 GeV heavier than the 



FIG. 2: The mass difference m^ — mfj 
as a function of m^, for the model hne 
described in the text. The soHd Hne is 
the model line with f^DM^^ = 0.11, and 
the shaded region denotes the approxi- 
mate region favored by the thermal relic 
abundance constraints. The model line 
is cut off on the left by the LEP2 Higgs 
mass constraint. 
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neutralino LSP. This is shown in figure [21 which plots the mass difference m^ — m^ as a function 
of nig for the model line. The bulge region where the stop-LSP mass difference is relatively large, 
with rrig between about 525 and 650 GeV, is characterized by having NiNi -^ ti due to ti exchange 
as the dominant annihilation effect in determining the dark matter thermal relic abundance. Note 
that varying tuq to obtain J7dm^^ anywhere within the allowed range 0.11 ±0.02 would not change 
the fact that m^ — tti^ < mw + rrn, for this model line. This stop-mediated annihilation region 
is continuously connected in parameter space to more fine-tuned models in which the ti, A^i mass 
difference is just right to allow efficient stop-neutralino co-annihilations, -f for rrig less than about 
525 GeV and greater than about 650 GeV. An important consequence of the larger stop-LSP mass 
difference in the dark matter annihilation-to-tops bulge region is that the LHC signal efficiency 
will be increased compared to the co-annihilation regions on either side, since the jets from the 
decay ti — > cNi tend to have higher pT- 

The superpartner production cross-sections are dominated by gluino and squark production. 
The next-to- leading order total cross-sections for this model line are shown in figure [3l computed 
using Prospino2 511]. The largest single source of supersymmetric events is pp — > tit^, which is 
of order 25 pb throughout the bulge region, and falls rather slowly with the gluino mass along 
the model line. However, this leads to the very difficult signal of two often low-py charm jets 
and little missing energy. I have checked that after realistic cuts to remove QCD and detector 
backgrounds (see for example [3, 153] )i the low efficiency of the iiil signal will lead to it being 
buried beneath the other squark and gluino sources, so it is not possible to infer the existence 
of the light top squark from this direct production process. The total gluino-squark associated 
production pp -^ gq plus gq*, summed over quark flavors, is of order tens of picobarns throughout 
the bulge region. Gluino pair production and (anti-)squark pair production both contribute of order 
10 pb in the bulge region, with the former falling somewhat more steeply with increasing mass. 
The production and decays of gluinos and squarks in this scenario should easily allow for early 
discovery lepton+jets+E'™''^'^ channels (see, for example, refs. [2, y, l53|] for comparable mSUGRA 
studies) at the LHC. 



* An even more fine-tuned stop-neutraUno co-annihilation region can also be found '5u\ in mSUGRA models. 



FIG. 3: The NLO production cross- 
section for superpartner pairs in pp colli- 
sions at y/s — 14 TeV for selected points 
along the model line described in the 
text, as a function of the gluino mass. 
Prospino2 51] was used. The most im- 
portant contributions, from stop pair 
production (iiii), gluino-squark pro- 
duction (gq and gq*), gluino pair pro- 
duction (gg), and squark pair produc- 
tion {qq and qq* and q*q*) are also 
shown separately. 
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Sleptons decouple from practical LHC physics in many compressed supersymmetry models, 
and in particular for the model line studied here. For example, for the model line point with 
Ml = 500 GeV shown in figure[Tl the total direct production cross-section of sleptons and sneutrinos 
before any cuts or efficiencies is only about 6 fb, compared to much larger backgrounds from WW 
production and other sources. As noted above, sleptons also extremely efficiently decouple from 
decay chains of heavier superpartners. Charginos and neutralinos (other than the LSP) do appear, 
but only in subdominant decay modes of the squarks. Their direct production rates are quite 
small compared to the gluino and squark rates. For example, for the model shown in figure [H one 
obtains a total of (129, 41, 32) fb for, respectively, {CfNj, CfCj, NiNj) production. Neutrahnos 
and charginos produced in association with gluinos and squarks adds another 500 fb. These rates 
are quite small compared to the 56 pb total gluino and squark production rate, and involve a wide 
variety of dissimilar final states without strong distinguishing features. Furthermore, these do not 
yield dilepton mass edges, as noted above. Unfortunately, finding out any information about the 
superpartners other than the squarks, gluino, and LSP from direct observation appears to be a 
daunting challenge at the LHC in this scenario. 

There are several ways of gaining information about the gluino and squark mass spectrum from 
the early discovery inclusive jets -|- leptons -|- S™'^^ signal, including for example mT2 and similar 
variables 5J, l55|] and the multiplicity of 6-tagged jets. However, the presence of non-negligible 
backgrounds that will have to be understood from LHC data puts these methods beyond the scope 
of the present paper. Instead, I will concentrate on tools that use the lower rate but potentially 
very low-background signal with same-sign dileptons. 



III. SIGNAL FROM SAME-SIGN LEPTONIC TOP DECAYS 



To define a signal for same-sign £^£'^bbjj + EJp^^^ events, I used MadCraph/MadEvent 56] 
for event generation, interfaced to Pythia 571] and then PGS4 J58t ] for detector simulation using 
CMS-like parameters. Events are selected by requiring the following from objects generated by 
PCS: 



FIG. 4: The number of LHC signal 
events with two same-sign leptons, two 
b tags, and two additional jets, per fb~^, 
after the cuts described in the text, for 
the model line described in section HIl 
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exactly two same-sign isolated leptons {i = e, fi) with pT > 20 GeV and |r/| < 2.4. 

at least two 6-tagged jets each with px > 50 GeV (with \r]\ < 1.75 required by PGS). 

at least two more jets with pT > 50,35 GeV (with \ri\ < 3.1 required by PGS). 

at least one pairing of each of the two leptons with a distinct ^-tagged jet, with each pair 
having invariant mass consistent with leptonic top decay: m{hi) < 160 GeV. 

^f?**^*^ > 100 GeV. 



(These cuts are very similar to those used in ref. 37[.) The uncorrected jet momenta from PGS are 
used. The reason for the cut on the bi invariant mass is that the parton-level kinematic endpoint 
is: 



m{be)r_ 



mt 



m 



ly 



(3.1) 



nominally about 153 GeV using mt = 172.7 GeV. I use the higher value of 160 GeV for the cut in 
order to partially take into account the effects of smearing of the 6-jet energies. The 6-tag is actually 
a heavy flavor tag, which in PGS has an efficiency for high-pT central jets of approximately 50% 
for true b jets, 13% for c jets, and 1% for g, u, d, s jets. Each of the leptons and jets are required to 
be isolated from each other by /S.R = y^{A(/))'^ + (A?])^ > 0.4. Also, muons that are not isolated 
from a jet are absorbed into the jet, if the summed px (excluding the muon itself) in a cone of 
AR = 0.4 around the muon exceeds 5 GeV, or if the ratio of the px in a 3 x 3 grid of calorimeter 
cells around the muon to the pT of the muon itself exceeds 0.1. 

The cross-section after these cuts for LHC collisions with ^/s = 14 TeV is shown in figure [H 
for points along the model line. Also shown are the two largest contributions, from gg production 
and gq production. The efficiency for the gg part of the signal is about 0.04%. Most of the gq 
contribution to the signal is due to production of a left-handed squark in association with a gluino, 
since right-handed squarks usually decay directly to the LSP and the corresponding quark. The 



cross-section after cuts is between 5 to 15 fb for gluino masses less than 640 GeV, corresponding 
to the bulge region where stop-mediated annihilation to top quarks dominates the dark matter 



annihilation in the early universe. For comparison, ref. [37| found backgrounds totaling less than 
0.5 fb, using very similar cuts (although a different event generation and detector simulation). 
Therefore, strong evidence for this source of supersymmetric events might be obtained with as 
little as a few fb~^, depending on the gluino mass, and to a lesser extent the squark masses. This 
of course presumes that the backgrounds can be well understood, and that wrong-sign assignments 
of lepton charges in e.g. ti production are indeed not large and irreducible. In the following, I will 
optimistically assume this to be the case, and neglect backgrounds. 

Note that the cross-section after cuts is actually lower for the lowest mass point in figure [3] with 
Ml = 425 GeV and rrig = 511 GeV than for the next-higher mass point with Mi = 450 GeV and 
rrig = 542 GeV. This occurs for two main reasons. First, the efficiency is lower for Mi = 425 GeV 
because of the much smaller stop-LSP mass difference, as noted in figure [2j Second, the Mi = 450 
point has much larger branching fractions for right-handed squarks to decay into gluinos, adding 
to the signal. 

[The contribution of right-handed squarks declines again for heavier masses, and is eliminated 
by kinematics for points on the model line with rrig larger than about 600 GeV. The branching 
ratios for dji -^ dg and ur -^ ug are (46%, 38%, 12%) and (21%, 16%, 5%), respectively, for the 
model line points with Mi = (450, 475, 500) GeV. For all other model line points shown in figure 
[3] and [31 the branching ratios of qr — > qg are negligible.] 

To give an idea of the characteristics of the signal events, figure [5] shows the E^^^^, lepton pT, 
and jet pT (with and without 6-tags) distributions for the events that pass the cuts, given 100 fb~^ 
of data for a representative model with Mi = 475 GeV {rrig = 569 GeV) . Clearly, raising the E^^^^ 
cut much farther above 100 GeV would have a significant unfortunate effect on the signal, even 
for heavier masses. The same is true for the subleading lepton and jet pr's- On the other hand, 
there is considerably more room to raise the cuts on the leading lepton and jet pt's without a 
huge effect on the signal, should that prove necessary to reduce backgrounds. Due to the practical 
difficulties that are anticipated in commissioning E^^^^ at the LHC, it is also tempting to consider 
dropping that cut altogether, since the same-sign dileptons and jet cuts alone might be enough to 
distinguish the signal from background. This may be, but figure [5] shows that the benefit accrued 
to the signal cross-section from relaxing the EJp^^^ cut below 100 GeV is limited, especially for the 
critical case of models with heavier gluinos, so for the purposes of the present analysis it will be 
kept. 

The frequency of heavy-flavor-tagged jets in the signal sample can help to confirm that the 
signal is really due to gluino pairs decaying to stops that in turn decay to charm quarks and 
LSPs. The number of events with 2, 3, 4 or 5 heavy-fiavor-tagged jets is shown in figure [6] for 
the point on the model line with Mi = 500 GeV, for 100 fb~^. Also included is the breakdown 
of these events into tight 6-tags as reported by PGS, with efficiencies for central high-pT jets of 
approximately 40% for true 6-jets, 9% for c-jets, and 0.1% for g,u,d,s jets. As a simple-minded 
check, one can assume that the m{bi) < 160 GeV requirement preselects only events that have 
the true 6-jets tagged, so that additional heavy-fiavor tags come from the true charm jets, and 
the numbers of events with 2, 3, 4, and 5 heavy flavor tags should be roughly in the proportion 
712 : n3 : 71-4 : ns = (1 — Pc)'^ : 2Pc{l — Pc) : P^ : 0, where Pc is the probability of a true charm 
jet to get a heavy flavor tag. Using n2 = 930 for the example in figure O and Pc = 0.13, one 
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FIG. 5: Representative transverse momentum distributions for i^f^bbjj + E^^^^ events from 100 fb~^ 
of superpartner production, after the cuts described in the text. The upper left panel shows the E™^^^ 
distribution (here without the E^^^^ cut) for two points on the model line with rrig — 569 and 702 GeV. The 
upper right panel shows the leading and subleading lepton px distributions for the model with irig = 569, 
and the lower panels show the leading and subleading 5-jet distributions (left) and the three leading jet 
distributions (right) for the same model point. 

would predict n^ = 278 and 714 = 21 and n^ = 0, in not unreasonable agreement for such a naive 
estimate with the actual finding of 71.3 = 299 and n^ = 37 and n^ = 1. This information can be 
used to clearly distinguish the scenario under study here from similar ones in which the stop-LSP 
mass difference is large enough to allow ti — > bWNi, or where strict minimal flavor violation leads 
to a dominant or competitive 4-body decay ti — > bff'Ni, either of which would lead instead to a 
parton-level same-sign dilepton signature of 



pp -^ l^e'Hbbbjjjj + Ef' 



(3.2) 



In this case, one would clearly expect many more events with 3, 4 and even 5 or more (since half 
of the hadronic W decays will result in a true charm jet) heavy-flavor tags relative to the number 
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FIG. 6: The number of signal events af- 
ter cuts in 100 fb~^ with exactly 2, 3, 
4 or 5 PGS heavy-flavor-tagged jets, for 
the point on the model line described in 
section|n]with Mi = 500 GeV, resulting 
in rrig = 596 GeV. The breakdown into 
numbers of tight b tags is also shown. 
The relative frequencies of heavy-flavor 
tags provides additional evidence for the 
g — *■ tti and ti — *■ cNi interpretation of 
the signal. 
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with 2 tags, compared to the situation in figure [6l Measuring the numbers of heavy-flavor tags 
within the i £ bbjj + i?™^^ signal sample can therefore establish whether m^ — m,^ < 85 GeV. 
Of course, the specific numbers for heavy-flavor tagging in the actual LHC detector environments 
might be quite different from those assumed here, but the principle should still apply. 

In the remainder of this paper, I will examine some strategies for obtaining information about 
the gluino, squark and LSP mass spectrum. Note that variables like 7tt.t2 SJ, |55| are hampered, for 
the same-sign dilepton event topology, by the inevitable presence of two neutrinos with unknown 
momenta in addition to the two LSPs in each event. I have therefore not attempted the difficult 
task of seeing whether this can give useful information when applied to the same-sign lepton sample. 
The definite absence of dilepton mass edges eliminates another commonly used tool [45|]-[49|] for 
reconstructing superpartner decay chains. Instead, I will consider mass estimators that use scalar- 
summed transverse momenta and single-lepton mass edges from visible gluino decay products. 



IV. MASS ESTIMATORS FROM SCALAR-SUMMED TRANSVERSE MOMENTA 

One of the most important efforts in a future LHC analysis of supersymmetry will be to obtain 
measurements, or at least estimates, of the superpartner masses. The purpose of this section is 
to consider observables that can serve as estimators of the masses of the superpartners produced, 
using scalar sums of the lepton and jet pt^s. There has been considerable effort in this area, often 
using the observables Ht and Mefj for events with jets and EJ^^^^. Here, I will study the prospects 
for using similar observables, but in the hopefully cleaner context of the £ £ bbjj + E^^^^ signal 
discussed in the previous section. To this end, consider four mass estimators defined by: 



Ma = Y.PTiJn)+ E PT{in)+Er^ 
n n=l,2 

Mb = Yl PTiJn) + Yl PTi^n) + Ef 
n=l, 2,3,4 n=l,2 



(4.1) 
(4.2) 
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Mc = J^PTOn) + Yl ^^(^")' (4-3) 

n n=l,2 

Mz? = J] PT(jn) + j; Pt(4), (4.4) 

n=l, 2,3,4 n=l,2 

where the jet labels are ordered by priji) > Pt{J2) > PriJs) > • • •• As usual, the idea is that the 
transverse momenta of the decay products should be approximately linear in the mass of the pair 
of heavy particles produced. The first observable, M^, simply sums over all visible object j?t's and 
the E}p^^^ . The second observable Mb is motivated by the ideas that the sum over only the leading 
4 jets should be less sensitive to theoretical uncertainties due to extra jets from the underlying 
event, and that the signal includes at least four quark partons. The other two observables Mc\d 
are the same as M^^b except that EJ^^^^ is not included. This is motivated by the fact that E^^^^ 
may be particularly difficult to obtain accurately, especially in early running of the LHC. 

Using the same event selection criterion as in the previous section, the distributions for these 
four mass estimators are shown in figure [71 for 100 fb~^ of data with three representative points 
on the model line defined in section [TI] with nig = 542, 596, and 675 GeV. Even from these coarse- 
binned distributions, it is apparent that the shapes of the distributions are distinguishable from 
each other. 

To determine a sharper empirical relation between these mass estimators and the superpartner 
mass scale, I performed an unbinned maximum likelihood fit to 100 fb~^ of generated events for 
each of 9 model points. Because the distributions of Mx (for X = A,B, C, D) are clearly far from 
Gaussian, better results are obtained by fitting them instead to the class of functions known as 
generalized inverse Gaussian distributions (with x = Mx)'- 



f{x) = -{x-xq) "exp 
n 



2{x - xq) J 



Here a, 6, and c are the fit parameters, and xq is the minimum of the distribution following simply 
from the jet and lepton px and £^™'^^ cuts. (In the present analysis, xq is equal to 325 GeV for 
X = A,B and 225 GeV for X = C, D.) The normalization condition 

f{x)dx = l (4.6) 

implies that 

ln(n) = ab + {1 - c)ln{a) +ln[2Kc-iiab)], (4.7) 

with Ki{z) the modified Bessel function of the second kind. The peak of each distribution, defined 
as the value where df/dx = 0, is then obtained as 



Mf ^'^ = XO + {Va%^ + c2 - c)/b, (4. 
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FIG. 7: Distributions for the mass estimators Ma,b.c.d defined by eqs. (|4.ip - (|4.4p . for 100 fb ^ of events 
for three models afong the model line described in section |TT1 with Mi — 450, 500, and 575 GeV. 

where a, h, c are set equal to their best-fit values ^^ after maximizing the log likelihood function. 

Performing a linear regression (with the variance for each model taken inversely proportional 
to the number of events found in 100 fb^^) gives relationships between the mass estimators M^'^ 
and the gluino mass: 



Mg = 1.693MP'^'' - 776 GeV, 
= 1.733MP''''' - 634 GeV, 
= 2.274MP'^'''' - 825 GeV, 



2.422MP^'''' 



676 GeV. 



(4.9) 
(4.10) 
(4.11) 
(4.12) 



* The best fits obtained in the following almost always turn out to have c very close to 3/2, corresponding to the 
special case known as an ordinary inverse Gaussian distribution (not to be confused with a normal Gaussian 
distribution). I do not know the explanation for this. 
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FIG. 8: The gluino mass is compared to the fitted peak values of the distributions of the mass estimators 
Ma,b.c,d defined by eqs. (|4.1[) - (|4.4p . Results for 100 fb~^ of events for each of 9 individual models along 
the model line arc shown as black dots, together with the best fit lines eqs. (|4.9p - (|4.12p . 

The comparison between the linear fits and the values obtained for the individual models are shown 
in figure [8j The smaller slopes of the Mg vs. M^^^ lines would seem to make them more useful as 
mass estimators than M^°^ , although this depends crucially on the presently unknown quality of 
the EJp^^^ determination. 

In general, the mass estimators might be expected to be roughly proportional to M^s = MsusY — 
mj^ /MsusY, where Msusy is a signal cross-section weighted average of the superpartner masses, 
in this case the gluino and left-handed squark masses. In models with a slightly larger rriq^ — rrig 
mass difference, the decay qn — > qg would be more important, and m,g^ would be weighted more 
strongly into MsusY- In the model line under study, the LSP mass and the gluino mass are very 
nearly proportional, and the squark masses are also tightly correlated with the gluino mass. In 
general, since the presence of the signal depends crucially on the Majorana nature of the gluino, 
this method should be useful to obtain a rough estimate of the gluino mass, albeit with some mild 
model assumptions. 



ENDPOINTS OF VISIBLE GLUINO DECAY PRODUCTS 



Another method that can be used to gain information about the superpartner masses is to look 
at the invariant mass distributions of identified visible products of the gluino decay. This method 
has already been extensively studied in ref. 37|, in a situation similar to the present one, but with 
a relatively much lighter top squark and LSP, and taking the squarks to be much heavier. In the 
model scenario under study here, the presence of a large component of gq — > ggq in the signal 
causes a significant additional source for confusion in identifying the jets following from the gluino 
decay. 

The parton-level kinematic endpoints from the decay g -^ tt\ followed by ti — > cNi and t -^ Mv 
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FIG. 9: The parton- level predictions for 
the distributions of the invariant masses 
m{Mc), m{bc), and m{£c), for the de- 
cay g -^ tti followed by ii -^ cNi and 
t -^ bW and W -^ iv. The masses 
used are nig = 596 GeV, rrif^ = 261 
GeV, and m^y = 201 GeV, correspond- 
ing to a point on the model line with 
Ml = 500 GeV. The endpoints of the 
distributions are at m{b£c) — 354 GeV 
and m{bc) = 283 GeV and m{ic) = 320 
GcV. 
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max + W.J — m,jy, 



m\ic' , ^2 ^2 



(5.1) 

(5.2) 
(5.3) 



where X{x, y, z) = x +y + z — 2xy — Ixz — lyz. Note that these endpoints are not independent; 
knowing any one of them yields the others, given the known masses of the top quark and W 
boson. (The widths of the particles, and the mass of the bottom quark, are neglected here.) The 
corresponding distributions [331 fo^ ^^^ model depicted in figure [D with mg = 596 GeV, m^^ = 260.5 



GeV, and rrijy = 200.8 GeV are shown in figure [9l Ref. 37| performed fits to the shapes of be and 



ic mass distributions, finding that the quality of the fits was made worse by leptons from taus in the 
top decays, among other effects. In the present case, there is a serious additional (and in practice, 
unknown) effect on the shape from wrong jet assignments due to the presence in some events of 
Ql ^ Qg, and to a lesser extent Qr ^ qg- Therefore, it is probably more robust to concentrate on 
the endpoints of the distributions. From figure [9] one sees that the m{ic) distribution is extremely 
shallow near the endpoint, making it very difficult to determine the endpoint from data. The m(bc) 
and m{bic) distributions are much steeper near their respective endpoints, so I will only consider 
them. It should be noted that different events contribute to the near-endpoint regions of these two 
distributions, even though the positions of the endpoints are algebraically related. 

To mitigate the problem of wrong jet assignments, I use a subset of events selected by the 
procedure described in section HIH with the additional constraint that the pairing of the two leptons 
with 6-tagged jets consistent with top decays [m{bi) < 160 GeV] is unique. (This reduces the signal 
efficiency by about a factor of 3.) For each bi pair, the putative charm jet is taken to be the one with 
the smallest value of m{bic) selected from among those with pT > 35 GeV. This selection means 
that far below the endpoints, there may well be many wrong assignments (both from extra jets in 
the underlying event and from jets produced in squark decays being assigned to the charm jet role). 



16 



but near the endpoints the assignments are made correctly with greater frequency. Results for 
the m{bic) and m{bc) distributions selected in this way are shown in figures [TOl and [TT] for several 
representative models. These distributions are seen to be roughly consistent with endpoints at 
the nominal positions, but wrong assignments and jet energy smearing leads to some events in a 
high-mass tail in each case. This can be seen to be particularly troublesome for the lowest-mass 
Ml = 425 GeV {rrig = 511) model point, where the small mass difference m^ — rrifj = 29 GeV 
means that the true charm jets often fail the px > 35 GeV cut.^ This exemplifies a more general 
difficulty. If the thermal relic abundance of neutralinos does not account for all of the dark matter, 
then the stop-LSP mass difference will be smaller than indicated in figure [21 for any given gluino 
masses. This can always lead to the problem of the charm jets having too small px and being 
replaced in the analysis by interlopers, leading to a distorted distribution and a tail above the true 
mass endpoint. To counteract this problem, one could use an independent check on the identity of 
the charm jet. In figures [TOl and [TTl the solid histograms show the portion of the signal for which 
the putative charm jet has a PGS heavy-flavor tag or contains a non-isolated muon (similar to the 
"soft muon" tag used in Fermilab Tevatron analyses). This information will clearly be more useful 
if the efficiency and purity of "charm tagging" can be improved. Although I will not attempt it here 
in the absence of a fully realistic detector simulation, one can imagine that a likelihood fit taking 
into account these effects could give measurements (or at least constraints) on these endpoints. 
However, figures [10] and [TT] show that precision may be difficult to achieve without either more 
data than 100 fb~^ or a better handle on charm jets. 

Unfortunately, even with such a measurement, most of these models are quite indistinguishable 
from each other using the endpoints or shapes of the distributions alone. This is because the 
endpoints are nearly independent of the mass scale defining the point along the model line studied 
here, as illustrated in figure [EJ In fact, for the entire range 545 GeV < rrig < 830 GeV, m{b£c)raa.x 
is within 10 GeV of 350 GeV for this model line. (For lower values of the gluino mass, the endpoint 
is lower, but its determination becomes much more problematic due to wrong jet assignments due 
to the smaller stop-LSP mass difference, as we have just seen.) It might at first seem surprising that 
the position of the endpoints does not scale with the gluino mass. The reason is that the scaling 
is counteracted by the factor of (1 — m^ /m? ) in the formulas eqs. (J5.ip - (j5.3p . which decreases as 
one moves to higher masses along the model line, because of the constraint on the stop-LSP mass 
difference coming from the dark matter abundance observation. 

Nevertheless, a successful determination of the endpoints will still be useful when combined 
with the information that the gluino decay signal is kinematically allowed at all. This is illustrated 
in figure [TS] which assumes that the m{bic) endpoint is found to be 350 GeV (or equivalently 
that the m{bc) endpoint is found at 279 GeV). The allowed line in the gluino mass vs. stop mass 
plane is shown, for various assumptions about the stop-LSP mass difference. For the signal to 
occur at all, one must have mi_^ — m^ < 85 GeV, otherwise the decay ti — > cNi would lose to 
the flavor-preserving three-body decay ti — > bWN\. (As noted at the end of section Hill this can 
be ruled out by counting the number of additional heavy-flavor tags in the events that pass the 
signal selections.) This means that for a given stop mass, the gluino mass must be above the solid 
line. The dashed lines show the gluino-stop mass relation for smaller values of the m^ — ?7ijy mass 
difference. [If the m{b(.c) endpoint is only constrained to be > 350 GeV, then the allowed regions 



^ I have checked that lowering this cut does not help significantly, because doing so also allows more interloper jets. 
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FIG. 10: The distributions of m{bic) for six points on the model hne described in section [IT] with Mi = 425, 
450, 475, 500, 525, and 550 GeV. The nominal endpoints, indicated by the dashed vertical lines, are given 
by eq. (|5.2p as, respectively, m(6^c)max = 273, 339, 350, 354, 350, and 345 GeV. The solid histograms show 
the portion of the signal for which the putative c jet has a heavy flavor tag or a soft muon tag. 



18 



Sioo 

o 



> 

CD 

CD 
o 

c 

CD 
> 
LU 



50 



"0 



o 
o 



100 
80 



> 


60 


CD 




CD 




o 

CM 


40 






en 








£_ 


20 


> 




LU 





o 
o 



60 



50 



.E 40 



M 



1 — ' — ' — ' — •- 
511 GeV 



gluino 

M^,„^ =196 GeV 

stop 

M^sp = 167 GeV 



X 




100 200 300 

Mass(bc) [GeV] 



3 n , 



400 



— 1 1 1 1 1 1 r- 



— I 1 1 1 1 1 1 1 1 — 

M , . =569 GeV 

gluino 

M , =257 GeV 

stop 

M^sp =190 GeV 




1 00 200 300 

Mass(bc) [GeV] 



400 



— I — I — I — I — I — I — I — I — I — I — 

622 GeV 



M 



gluino ' 

M , =263 GeV 

stop 

M^gp =212 GeV 




100 200 300 

Mass(bc) [GeV] 



400 



Sioo 

o 



> 

CD 
O 

o 
£i 
To 

CD 
> 

LU 



50 



f 



M , . =542 GeV 

gluino 



stop 
LSP 



246 GeV 
179 GeV 




80 



o 
°60 



50 



o 




o 






4U 


_c 




> 30 


CD 




CD 




o 


20 






CO 








0) 


10 


> 




LU 





100 200 300 

Mass(bc) [GeV] 



400 



— 1 1 1 1 r- 



~r 



Mgiuino = 596 GeV 
Mstop = 261 GeV 
M^sp = 201 GeV 




100 200 300 

Mass(bc) [GeV] 



400 



— 1 1 1 1 r- 



Mgiuino = 649 GeV 
Mstop = 267 GeV 
M^sp = 222 GeV 




100 200 300 
Mass(bc) [GeV] 



400 



FIG. 11: The distributions of m{bc) for six points on the model line described in section |TT] with Mx = 425, 
450, 475, 500, 525, and 550 GeV. The nominal endpoints, indicated by the dashed vertical lines, are given 
by eq. jO]) as, respectively, m(6c)max = 200, 268, 279, 283, 278, and 274 GeV. The soHd histograms show 
the portion of the signal for which the putative c jet has a heavy flavor tag or a soft muon tag. 
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FIG. 12: The parton-level prediction for 
the endpoints of the m{Mc) and m(bc) 
distributions, for points along the model 
line described in section |TT1 as a function 
of the gluino mass. 



FIG. 13: The relation between the 
gluino (g) and the lighter stop (ii) 
masses following from eq. ()5.3|) . for 
the case that the high endpoint of the 
m{bic) distribution is taken to be 350 
GeV. The different lines correspond to 
TOj^ - m^^ = 30, 45, 60, and 85 GeV. 
The last is the maximum mass differ- 
ence that allows ti — > cNi to dominate, 
as required by the signal. 
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are above the indicated lines.] Now, combining this information with an estimate or upper bound 
on the gluino mass from the production cross-section or from the observables of the type Ma,b,c,d 
described above would allow a determination of ranges in which the gluino, stop, and LSP masses 
must be. 



VI. OUTLOOK 



Compressed supersymmetry with top-squark mediation of neutralino annihilation in the early 
universe presents both challenges and opportunities for the LHC. Although early discovery should 
not be a problem because of the low mass scale, the sleptons and the charginos and neutralinos 
(other than the LSP) may very nearly decouple. In the model line studied here, for example, it is 
very difficult and perhaps impossible for the LHC to be able to say anything about them. With 
sufficient integrated luminosity, one may be able to discover stoponium through its diphoton decays 
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39l . |59|], giving a uniquely precise measurement of the top-squark mass. This would provide an 
important absolute reference point for determination of the other superpartner masses. However, 
it requires that the top-squark mass is not too large. 

Other than stoponium, the most distinctive signature may be the same-sign leptonic top-quark 
decays that come from gluinos (or squarks decaying to gluinos). In this paper, I have studied the 
prospects for learning about the gluino, top-squark and LSP masses from these events. 

First, the scenario is distinguishable from similar ones with a larger stop-LSP mass difference 
by using the frequency of additional heavy-flavor tags in the i^i^bbjj + EJ^^^^ events after cuts. 

Observables obtained by summing over scalar pt's and E^^^^ within this relatively clean sample 
will provide estimates of the effective superpartner mass scale, which is always strongly correlated 
with the gluino mass. This can be compared with the estimate obtained from Ht and M^q 
distributions in the larger inclusive jets-l-i?™ ^'^'^ sample. 

The determination of invariant mass endpoints is somewhat more problematic, due to the per- 
nicious effects of interloping jets (both from squark decays and from the underlying event) being 
confused with the charm jet in the analysis. In the actual LHC analysis, this can probably be 
enhanced by using heavy-flavor likelihoods on an event-by-event basis to help choose the correct 
charm jets. These endpoints do not provide unambiguous information about the superpartner 
masses, even within the confines of the single model line studied here. However, when combined 
with the information that the decay ti — > cNi dominates, this information can also be useful to 
constrain the model. Clearly, heavy flavor tagging will be crucial in this effort. Also, if one is will- 
ing to assume that the thermal relic abundance of dark matter is due entirely to neutralino LSPs 
without fine-tuning, then the resulting stop-LSP mass difference should be large enough to more 
sharply define the endpoints. Conversely, a confirmation of this scenario would help to establish 
the supersymmetric interpretation of the dark matter. 
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